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Abstract 

A data set comprising 110 spreading rates, 78 transform fault azi- 
muths and 142 earthquake slip vectors has been inverted to yield a new 
Instantaneous plate motion model, designated RM2. Tlie model represents 
a considerable improvement over our previous estimate KMl (Minster, 

Jordan, Molnar and Haines, 1974) . The mean averaging interval for the 
relative motion data has been reduced to less than 3 My. A detailed 
comparison of EM2 with angular velocity vectors which best 
fit the data along individual plate boundaries indicates that RM2 per- 
forms close to optimally in most regions, with several notable exceptions. 
The model systematically misfits data along the India-Antarctica and 
Pacific-India .plate boundaries. We ^hypothesize that these discrepancies 
are manifestations of internal deformation within the Indian plate; the 
data are compatible with NW-SE compression across the Ninetyeast Ridge 
at a rate of about 1 cm/yr. RM2 also fails to satisfy the EW-trending 
transform fault azimuths observed in the FAMOUS area, which is shown to 
be a consequence of closure contraints about the Azores triple junction. 
Slow movement between North and South America is required by the data set, 
although the angular velocity vector describing this motion remains poorly 
constrained. The existence of a Bering plate, postulated in our previous 
study, is not necessary if we accept the proposal of Engdahl and others 
that the Aleutian slip vector data are biased by slab effects. Absolute 
motion models are derived from several kinematlcal hypotheses and compared 
with the data from hotspot traces younger than 10 My, Although some of 
the models are inconsistent with the Wilson-Morgan hypothesis , the overall 
resolving power of the hotspot data is poor, and the directions of absolute 
motion for the several slower-moving plates are not usefully constrained. 
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’ Introduction 

Present-day plate motions can be modelled using systematic 
Inversion methods. In our initial study (Minst ar et al. , 1974, 
referred to as Paper I), a linearized least-squares algorithm was 
formulated and applied to an extensive, globally distributed data 
set. Angular velocity vectors for eleven ‘major plattss were estimated 
from these data, and this model was designated Relative Motion 1 (RHl) . 
The Caribbean plate was subsequently added to this model by Jordan (1975). 
Revisions and additions to the data set were begun in 1975, and an 
Interim model was derived (Jordan, Minster and l-Iblnar, 1976). 

We present in this paper a new relative motion model, EM2, based 

on a much improved data set. Consistent irith our previous work, 

' * 

we have attempted to obtain a simple model compatible with the available 
high-quality observations of relative motions. Only relative motion - 
data whiph involve at least one oceanic plate have been used, since the 
data from intracontinental environments exhibit complexities not easily 
described in terms of rigid plate kinematics (e.g. Molnar and Tapponnier, 
1975). We have not attempted to model the complex tectonics 
of the western Pacific (e.g., the Philippine plate), because little 
klnematical information is available concerning behind-the-arc spreading, 
and the assumptions fundamental to a simple plate model (e.g. triple- 
junction closure) may not apply. 

The value of any model can be judged by its predictive capability 
and by its ability to withstand the test of new observations. In this 
respect the success of our original model RMl has been mixed. For 
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example, the relative motion between the Worth Ameriran and South 
American plates was predicted by EKL entirely on the basis of data 

■V 

from other plate boundaries. Although no data yet exist which confirm 
directly the existence of such relative motion, the model implies that 
a component of NS convergence exists between the South American and 
Caribbean plates (Jordan, 1975) . It appears that some convergeftce is 
indeed required by recent studies (Talwanl et al. , 1976; Rial, 1978). 

On the other hand, RMl failed to satisfy an extensive set of new 
data collected in the South Atlantic Ocean (Forsyth, 1975; Sclater 
et al. , 1976) . The investigation of this failure is an important 
aspect of this study. We show that RMl incorrectly predicts the 

plate kinematics in the South Atlantic because the presently available 

* 

data are inconsistent with the plate geometry assumed in deriving EMI. 

We demonstrate that this inconsistency can be remedied by postulating 
the existence of internal deformation' with the Indian plate, although 
alternate explanations are possible. 

Other problems with the EMI model have beeti noted. (Jordan et al. , 1976) 
The well-mapped fracture zones in the FAMOUS area yield an apparent 
azimuth for Africa-North America motion that is due east (Macdonald and 
luyendyk, 1977), whereas EMI predicts an azimuth of S79E, parallel to 
the general trends of the nearby major transform faults (e.g, the 
Oceanographer T. F.). 

In EMI the slip vector data from the North Pacific were modelled 
-using a Bering plate whose motion differs from that of North America, 

Engdahl et al, (1977) have demonstrated that the focal mechanisms from 
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this region can be affected by slab structure, perhaps biasing the 
observations. They have suggested that corrections for this bias may 
eliminate the need for a Bering plate. 

These and other problems are examined in this paper. 





The Revised Datfa Set 

The 330 data used In this study are listed in Table 1. The data 
locations are shown in Figure 1, delineating the major plate boundaries. 
These relative motion data comprise 110 rates of sea floor spreading 
derived from magnetic anomaly profiles, 78 transform fault azimuths 
and 142 earthquake slip vectors. In compiling and editing this data 
set, we have generally followed the guidelines in Paper I. In particular 
we have excluded data from diffuse plate boundaries, specifically 
continent-continent boundaries. Therefore the details of Asian and 
Indonesian tectonics are not represented by our model. 

Rate data have been determined directly from published magnetic 
anomaly profiles using the time scale of Talwani et al. (1971) • In 
Paper I, anomalies 3 and 5 were generally used to estimate rates; we 
thus averaged the plate speeds over the last 5-10 My. In this study, 
we have redetermined the spreading rates using anomalies 2 and 2’ in 
every instance, except for a few slow-spreading profiles where the 
anomalies out to 3 were employed. Hence, the mean averaging interval 
for the rate data is less than 3 My. In most cases the rates were 
determined by comparing the corrected profiles with synthetics, 
generally those published by the authors of the 'original observational 
study. However, for the anomaly profiles along the Pacific-Antarctic 
Ridge (Molnar et al. , 1975), we generated our own synthetics. For tlie 
several studies where a direct Inversion for magnetization was made 
(Hacdonaldt 1977; Macdonald and Holcombe, 1978; McGregor et al. , 

1977), the original authors’ results were used directly. 




In Paper I, the directions of plate motion implied by earthquake 
focal mechanismE were estim^tfd by projecting the. slip vectors piitp a 
horizontal plane. Although this procedure is almost universally, adopted 
- it is only approximately copr.ect for shallow thrust events, in subduction 
..zones with obl.ique convergence, and it can introduce a sligh.t bias. 

.In this study, the more exact procedure of ‘rotating the slip vectors 
into the horizontal plane was employed for earthquakes along inclined 
seismic zones. This problem is discussed Ln the Appendix. 

The. m.ost, precise estimates of. relatiya motion direction are the 

azimuths of well-mapped transform faults. In determining these azimuths 
we have used detailed bathymetric surveys where available,., relying on 
contours which cross charted ship tracks. .Interpretive diagrams have 
been avoided to minimize the feedback between -data and plate tectonic 
. models. . ■ - 

The .uncertainties listed in Table 1 are based on a case by case 
subjective evaluation of the data quality. They are used to weight the 
..data in the inversion algorithm and to derive estimates of the. uncer- 
. tainties in the model. parameters. Although we have attempted, to use a 
consistent set of . criteria. in assigning these errors, the estimates are- 
nevertheless crude indicators of data quality. With this in mind, we 
have adopted a conservative .stand and have deliberately overestimated 
^ these uncertainties. This bias is apparent in Figure 3, where it is 
. seen that the sample standard deviation of the normalized residual 
distribution is significantly less than its expected value of 1. 
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Model 8H2; General Pescrlptioa 

Inversion of the data was performed using the linearized, iterative, 
weighted least-squares algorithm described in Paper X. Our extensive 
experience in applying this algorithm to the plate motion problem has 
demonstrated to us its effectiveness. Although the algorithm involves 
the linearization of a non-linear problem, convergence has always been 
rapid and no difficulties associated with local minima have been evident. 
The uncertainties in the model parameters derived from the linear theory 
have proven to be effective measures of the errors induced on the model 
by errors in the data. 

The inversion algorithm has been applied to the data set listed 
in Table 1 to obtain an 11 plate model, designated Kfl2. The plate 
geometry is identical to that of RMl, except that the Bering plate 
has been recombined with the North American plate and a Caribbean plate 
has been added, BMl, supplemented with the CARB-NOAM angular velocity 
vector derived by Jordan. , (1975) , was used as a starting model in the 
inversion algorithm. Convergence was attained in five iterations. 

Model RM2 is specified in Table 2 by its geohedron (McKenzie and 
Parker, 1974), Although a more compact specification is possible, this 
format conveniently provides an explicit relative rotation vector for 
each plate boundary. The Rtl2 geohedron is illustrated in Figure 2. 

In the notation of Paper X the quantity minimized by the fitting 
■procedure is the variable 

. I; 

* i=l „ 2 
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where N »* 330 is the total number of data. The elevt‘ii plate model is 
specified by 30 parameters. If the data were normally distributed and the variances 
Were exactly known, would be chi-square distributed with 300 degrees of freedom, 
and a sample value would lie in the interval (300 ± 49) 95% of the time. 

The value of for ^2 is 109, almost a factor of three less than its 
expected value. Thus, the data are fitted significantly better than they 
should be if their assigned uncertainties were corref^t. 

This fact is also evident from the histograms of normalized 
residuals plotted in Figure 3. The sample variances of these distribu- 
tions are about 1/3 their expected value of unity. This discrepancy 
could be corrected by uniformly reducing the standard errors assigned 
to the data by a factor of 1//3. Such a reduction would not change the 
model but would decrease the derived model uncertainties by the same 
factor. However, to be conservative we have retained the larger 
estimates of uncertainty. 

It can be seen from Figure 3 that the distribution of normalized 
residuals for the slip vector data departs from the assumed gaussian 
behavior in another manner; the distribution is skewed towards negative 

w 

values. Much of this skewness is attributable to the predominantly 
negative residuals exhibited by the slip vectors from the Aleutians 
and the Kurils, a feature discussed in more detail below. 

Because the data set is large and because the geometry of the 
problem is complex, the performance of EM2 cannot be fully described by 
these simple statistics. A complete assessment of RM2’s success in 
explaining the observations requires that each data. subset pertaining 



to stn individual plate boundary be considered separately. For a large 
number of plate pairs, a relative rotation vector, oi* at least a “best- 
fitting pole" (BFP) , can be determined from that data subset alone. 

These vectors and poles have been obtained by inversion and are listed 

V 

In Table 3. The corresponding BFP's are shown with the RMl and RM2 
poles on Figures 4-6. The differences between these poles and those 
for KM2 measure the constraints imposed on RM2 by the simultaneous 
inversion scheme. These differences are not large, which is evidence 
that MI2 performs close to optimally in most regions. Notable exceptions 
involve the INDI-ANTA, INDI-PCFC and AFRC-NOAM poles, discussed below. 

The estimated model uncertainties O., O are much smaller in 

8 <p w 

Table 2 than in Table 3. This is, of course, a direct consequence of 
tie self-consistency constraints inherent to the rigid plate model, 
as discussed in Paper I. An impressive example of this behavior is 
provided by‘ the COCO-FCFC rotation vector, which is heavily constrained 
by two triple junction closure conditions; these constraints reduce the 
nominal uncertainty of the rotation rate by a factor of four. 

It should be emphasized that the uncertainties in the model para- 
meters given in Table 2 correspond to marginal distributions. A 
complete description of the mode] uncertainties | including the various 
error cross-correlations, requires the specification of a 30 ^ 30 
(symmetric) variance matrix. A more complete discussion of this point 
Is given in Paper I. 

Listed in Table 1 are quantities which we have termed "data 
importhnees." As defined in Paper I, they are the diagonal elements 



of ^11 orthogonjil projection operator tlie_ datji spa<:e, ^nd_are_. 
indicative of th^ distribution pf_ Informat ion. among the (Pappr I; 

Mn^t^er et aL. , 1977^ Importances are additive and sum to .the number 
of inverted parameters, 30 in the. case of KM2. They depend on the 
geometry of the data set, and- on the data .uncertainties, but not on 
the actual values of the data. The final model depends heavily on the 
most Important data and is robust with respect to thu least important 
data. • 

Cumulative importances for individual plate boundaries are listed 
by data tjrpes in Table 2 for EM2 and in Table 3 for the best fitting 
vectors. The cumulative importance for all slip vector data is only 
4.6, compared with 11.1 for the transform fault azimitths, despite the 
fact that the former outnumber the latter by nearly 2:1. This reflects 
the lower uncertainties — by a factor of two to three — generally assigned 
to transform fault data. The most important datum (0.95) is the rate 
across the Mid-Cayman Rise (Macdonald and Holcombe, 1978); alone, it es- 
sentially determines the relative speed of MOAM-CARB. l^hen the entire data 
set is considered, 50% of the cumulative importance is associated with 
the 49 most important data, and only 10% with the 151 least important 
data. Importances are very useful for a detailed comparison of data 
and models, as illustrated in the next sections. 


This discussion is devoted to a detailed evaluation of BU2 on a 
region-by~region‘ basis . The fit of RMl and EM2 to the data for 
individual plate boundaries is illustrated in Figures 7-20. The data 
and model values are depicted as residuals with respect to the best- 
fitting angular velocity vectors and poles listed in Table 3. .Baselines 
provided by the best-fitting vectors remove the large variations in the 
data functionals due to geometrical compliaclties and allow the models to 
be plotted as smooth lines on the diagram.s. More importantly, the 
deviations from the locally best-fitting parameters required by closure 
conditions are readily apparent. 

The Pacific-Worth America Boundary. It was concluded in P.aper I . 
that the slip-vector data along the Aleutian-Kuril trench system are 
not consistent with the WOAM-PCFO relative motion inferred from data 
in the Gulf of California and in the northwest Pacific. We suggested 
that this inconsistency was diagnostic of deformation of the North 
American plate, and attempted to model It by including a hypothetical 
Bering plate in k£i. However, the BiiRT-PCFC pole was determined by only 
ten slip vectors. Engdahl et al. (lf>77) pointed out that our data were 
a poor representation of the earthquake population along the trench and 
that the slip vector orientations for individual events in the vicinity 
of 175**E could be significantly biased by the laterally heterogeneous 
seismic velocity structure of the downgoing slab. In the present study 
the number of data along this trench system has been Increased to 27, 
Including 15 high quality slip vectors from the Kuril-Kamchatka Arc 








recently published by Btamter and Mualchin CI976) * because of Che evi- 
dence for bias duerco slab structure presented by Enjtdahl et al. (1977), 
we assig'ne'd large uncertainties (±20°) to the data lying between 165°E 
■ ■ "and 165°W longitude. It can be seen from Figure 7 that these data are 
in fact systematically misfit by EM2 and the BFP in the direction ob- 
served in Paper I and predicted by the model of Engdabl et al. (1977), 

"On the other hand data from the Kuril -Kamchatka Ape are fitted by the 
model without difficulty, consistent with the conclusion of Engdahl et al. 
(1977) that slip vectors in this region are not likely to be significantly 
biased by slab structure. Since the fit of the data elsewhere along the 
boundary is satisfactory (Figure 7), we conclude that there Is little 
evidence for deformation within the Korth American plate of the sort 
hypothesized in Paper I. . , 

The East Pacific Rise. The data set for the COCO-PCFC boundary 
includes a redetermination of the Siqueiros T.F, azimuth from revised 
bathymetry (Rosendahl, 1976). R>!2 performs very well along this 

boundary and constitutes a slight improvement over KJtl (Figure 8) , 

The data set for the NAZC-PCFC boundary has been significantly 
revised and augmented, especially the rate data set. Between 6°S and 
12°S, the magnetics are poor and the data relatively scattered (Figure 9) , 
as might be expected for. east-west profiles In the vicinity of the mag- 
netic equator. Nevertheless, Rea's (1976t>, b) data indicate a lower 
rate than used in Paper I. Herron’s (1972) profile at 19°S is easily 
readablv“, despite the small size of the published figure, but the 
■ bathymetry indicates that a fracture zone may be crossed to the west of 
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the ridge. Thus, the western part of the profile is suspect "beyond 
anomaly 2, and we assigned a large uncertainty to the measurement. A 
sequence of high quality profiles at 20°S has been discussed by Rea 
and Blakely (1975) . Since their published profiles are rate adjusted 
and could not be remeasured, we adopted their estimated spreading rate 
(16.1 cm/yr) and assigned it an uncertainty of 0.6 cm/yr, a conservative 
value in view of the datum’s quality. However, this rate is less than 
that obtained at 19°S and is not fitted well by the model. It is also 
difficult to reconcile this rate with the comparable rates much further 
north and a higher rate to the south: the profile at 28°S (Herron, 1972) 

yields a rate which exceeds 17 cm/yr. 

The azimuths along the NAZC-PCFC boundary have been much improved 

r • 

by the recent bathymetric studies of Mammerickx et al. (1975) and 
Lonsdale (1977, 1978). However, the position of the NAZC-PCFC pole 
has not been significantly altered by these revisions; the BMl and RH2 
poles, and the BFP, lie very close together, well within the RM2 error 
ellipse. . 

The Galapagos Spreading Center. The rate data along the COCO-NAZC 

boundary are taken from the study by Hey (1974) . We also included a 
.• * 

good deep-tow profile published by K1 itgord and Mudle (1974) . As 
seen in Figure 10 and in Table 1, the data along this boundary are 
internally consistent. A particularly satisfying feature is that the 
recent bathymetry of Lonsdale and Klitgord (1978) clearly requires the 
COCO-NAZC pole to lie north of the equator; the transforms at 84.5°W 
and 85.3°W trend east of north. The implied shift from the RMl pole 


position is in complete concordance wi^ the shift dictated by the lower 
spreading rates along the NAZC-PCFC boundary. It should be noted that 
the strike of the Panama T.F. is very consistent with this new pole 
position, a point we shall return to in the nest section. 

The Chile Rise. The slower opening rate along the NAZC-PCFC 
boundary also affects the motion along the Chile Rise. In particular, 
the EM2 rate is considerably less than the 7.6 cm/yr estimate derived 
from the profile of Klitgord et al. (1973) , which we consider to be the 
best rate observation along this boundary and is the only value included 
in the data set. However, the RM2 rate if- between this value and the 
lower estimates of Morgan et al. (1968) and Herron and Hayes (1969) . 

Eastern Pacific Subduction Zones. Strongly coupled to the opening 


of the East Pacific Rise are the convergence rates and directions along 
the Middle American and South American trenches. We have adopted a 
set of slip vectors estimated by Stauder (1973, 1975) and Abe (1972) to 
represent the direction of subduction in South America, The residuals 
for these data show a slightly negative trend, although Abe's (1972) 
well determined solution has a large positive residual. The negative 
trend could be eliminated by increasing the rate along the NAZC-PCFC boundary. 
However, the COCO-NOAM and COCO-CARB slip vectors also exhibit this 
negative residual trend, and the possibility that these data are biased, 
like the Aleutian slip vectors, cannot be discounted. In any case, the 
scatter in the data is large, the average misfit is small and the data 


Importances are low; hence, any bias will not sigciificantly affect the 







- sThedBaclElc-rAtitraretlc ■Rt^dg^/ ~=ffi~cau5iP£t ~£^ s^tctiystone of the 
globaimodat, particular :atrtcntioii’vas d6yofBd“to the PCTC-ANTA boundary. 

The data -along this boundary are of 'sufficient -Tvuinber ‘ and quality to 
provide isignif leant 'coupling, via the Antarctic plate, among the plates 
in the Pacific and the plates tilth boundaries in the South Atlantic and 
Indian Oceans, The configuration of the PCFC-AHTA boundary has been 
Investigated by Molnar et al. (1975) , and our data set is based primarily 
on this study. Since these authors did not use synthetic magnetic . 
profiles, we computed synthetics and reinteirpreted the magnetics, A 
significant component of apparently asymmetric spreading is observed 
on many profiles (Molnar et al, , 1975; Stein et al. , 1977), so the 
rates were estimated only from pairs of corresponding anomalies on 
both sides of the axis All meanurements were based on anomaly 
or younger anomalies. Transform fault azimuths were derived from the 
bathymetry ' of Molnar et al, (1975), but estimates were obtained from 
ship track crossings rather than their interpretive map. It is clear 
from Figure 11 that RM2 is very close to the best-fitting vector and 
represents an imprjDvemont over RMl in this region. The difference in 
the RMl and KH2 poles is mainly attributable to the southwesternmost transform 
fault, an important datum (^^= t>.25) not included in Paper I. Some internal 
inconsistency of unknown origin is evident in the rate data (Figure 11) : 
the rates are greater in the middle of the boundary than those required 
by the rates at the ends of the boundary, nevertheless, most of the 
data are fitted within their uncertainties, and the relative rotation 
vector is one of the best determined in the BM2 geohedron. 
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'^r:Thfe".lhdia-Paelflc boundary*! The data used alang 'this boundary, 
*=i^censlstlng entirely: of -earthquake sllp'^eetots, “are'She* same as-:±n=“ 
-r.:paper I, but-the data :nGrth:of"25°S:tter6:eliininated: because of documented 
“~“b6hind-the~arc' spreading in‘ the"Lau basin' (e.g. : Lawyer •• et~ al , 1975) . 

: Nevertheless, the geometry is such tliat-a BFP could be determined -from 

!- :the 14- remaining slip vectors (Table 3) i We observe ' 

that this best-fitting jiole is almost identical to the pole determined 
by Falconer (1973) exclusively from seismicity data ?long the Macquarie 
: Ridge, a completely independent data-set. However, as seen in Figures 

*--4-and 12, -both RMl and RM2 differ significantly from this pole, a direct 
result of -requiring closure around tlie INDIrPCFC-ANTA triple 
•Junction. Consequently, the global models are a poor fit to the 
southernmost slip vectors, determined by Banghar and Sykes (1^69) . 
Furthermore, these models predict a significant component of compression 
across the Macquarie Ridge system, in disagreement with the hypothesis 
: .-of Falconer (1973) that this segment is a strike-slip fault. We 
-••strongly suspect that -these inconsistencies result from internal 
deformation -within the Indian plate (sec below) , ! . . .... 


Motions about the Azores triple ,j unction. The plate boundaries 
which form the Azores triple junction are individually well constrained. 
Figure 13 is a residual plot for the northern Mid-Atlantic Ridge data. 
The longitude of the EURA-NOAM pole is reasonably well fixed by the 
precise azimuth data along the Charlie-Glbbs T.F. and a number of fault 
plane solutions in the Arctic, but its latitude is more uncertain. Both 


the RMl ^ole ' and near th^ mouth of -the JLetia 1M.ver, the _ 

position ^ost 'compatible with-tht;, rate data.- -The BM2 pole_is several 
degrees furiher'south" (65;8°N, -132.4°E),-imd^lts_fit to.the rate dat^ 
sciith of~60^N is hot as good. -However, tMs pole more consistent, 
with the conclusions reached by .Chapman, and Solomon (1976} injfcheir 
study of northeast Asian tectonics. 

The -data set along the Azorcs-GSbraltar Line is considerably 
improved "over our previous study. We deleted the datum east of 
Gibraltar, because of its probable involvement with the Alboran plate 
(Andrietix et al. , 1971) , but added three new. slip vectors west of 
Gibraltar, The most important addition,^ however,. is the asimuth of the 
Gloria T.F, (^= 0.783), well defined by Laughton et al, (1972) and Laughton 
et al. (1975). This datum places *a strong constraint on the longitude of the 
AFRC-EURA pole. Although the individual slip vectors are not particularly well 
determined, their variation from NW compression on the east to SW extension on the 
west requires that the pole be not fur south of the boundary, a conclusion 
established by McKenzie (1972). As a result, the pole is very tightly 
' tonstrained, and the RM2 solution is very close to the BFP (Figures 6 
and 14) , 

-• The data set south of the Azores on the Mid-Atlantic Ridge has also 
been improved. Several special studies' have yielded much better mag- 
netics, and these imply a significantly lower rate during the last 3 My 
..than used in Paper I. The azimuth data along the AFRC-NOAM boundary 
have also been revised. In Paper I, the general trends of the Oceanographer 
T.F* (S77E) and the Atlantis T.F. (S81E) were used and were well fitted 


by EMX. In tbe present data set, these aslmuths have been deleted and 
replaced by the azimuths of transforms A(S88E) and B(S89E) in the FAMOUS 
area (Macdonald and Luyendyk, 1977). The difference between the azimuths 
of the major transform faults and transforms A and B has been attributed 
to a change in the direction of plate motion within Che last 5 My 
(Macdonald, 1977; Fox et al. , 1978; Atwater and Macdonald, 1977). 

A slip vector showing east-west motion on the Oceanographer T-F. 

(UdfCas et al. , 1976), supporting this coniilusion, has also been included. 

The revised data along the AFRC-NOAM boundary are internally 
consistent, as indicated by the performance of the best-fitting angular 
velocity vector, but the AFRC-NOAM azimuth data are poorly fitted by 
RM2 (Figure 15). It is clear that the misfit is forced by the closure 
condition about the Azores trip'le junction. To satisfy the 'triple 
junction condition, the AFRC-NOAM pole must be on the great circle 
connecting the EURA-NOAM and AFRC-EUEA poles (Figure 4 & 6) , The BFP 
Is not; it lies to the west near the northeastern tip of Greenland, 
as required by the revised azimuth data. The triple junction great 
circle cannot he shifted to include the AFRC-NOAM BFP without completely 
misfitting the data along one or both of the other boundaries. For 
example, any good fit to both the AFRC-NOAM and EURA-NOAM data sets 
yields an. AFRC-EURA pole that is much to the west of the RM2 pole and 
implies compressive motion along the entire Azores-Gibraltar Line, a 
prediction in flagrant disagreement with the observed earthquake mechanisms , 
Hence, the RM2 solution is significantly different from the AFRC-NOAM BFP. 
ilie.RMl and RM2 poles are each included within the other* s 95% confidence 
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elllp^a. dj.yeciJLp.iis^ o£- -ABR G-NQAM. mot^fm r r-' 

iaat_ch ±he^ pbseCTTjid A®n.eia.l .the. .Oi:fianograpl^ but whlcji : 

Elis fit th,e_ azimuths. . of ttausformf^ A-and B-Ly about- 

A pos.sible .explanation for this discrepancy concerns the iway- the 
RMZ, data set averages over time. It is conceivable that the east-west 
trends observed in the FAMOUS regionr are so "recent that the pole shifts 
required by this reorientation are not represented in. the data from the 
other, plate boundaries. 

Howeverj we believe tha.t this explanation can be rejected. The 
location. of the great circle connecting the EURA-HOAM and APRC-EUEA 
poles is fixed by truly ’’instantaneous”, dataj i.e., the slip vectors 
in the North Atlantic and along the AzoresrGibraltar Line. Therefore, 
the conflict is among data which involve little or no time averaging. 

Perhaps the east-west transforms observed in the FAMOUS area are 
not unbiased indicators of AFRC-KOAM motion.*^ This would be the case, 
for example, if these short, fault segments were ’’leaky” in the sense of 
Menard and Atwater (1969); i.e.,if a component of extension existed 
across these faults. For this explanation to be correct, the rate of 
opening normal to the faults would, have to he about O.A cm/yr. Although ' 
the field data do not appear to support this hypothesis (Detrlck et al. , 
1973; ARCYANA, 1975; Ghoukroune et al, , 1977), the ability of these 
studies (as well as ours) to resolve such a component is an open question. 

The incompatibility of the .FAMOUS trends with the RM2 model remains 
problematic. It is Interesting to note, however, that the EtI2-predicted 
.azimuths are essentially perpendicular to the rise-crest segments in the 


FAMOUS area. 



,r.* : 1 c Amaglcast ~ one - plate pr* two? A-' Ttajb^ conclusion^ of*-Pape r I 
:iBfas-that-slgiiiflcantrrelative motion existrs.between 'Worth. and South 
;Ajnerica,c-The*presentTstudy~supports this conclusion, although direct 
observationalrevidencB' for* NOAttrSOAM motion is still lacking. An 
inversion of the global data set Was perf(»nned with :he Americas grouped 
into a single plate,~ This model was rejected because • it does not satisfy 

r 

the relative motion data in the AtlantlCi.. In particilar:.’ 

(1) The rate data along thei AFRC-NO/\M boundary are misfit, model 

r values. being 0,4 cm/yr too low. • / ' 

t 

(2) The azimuths along the APRt:-S0AJl boundary vield systematically 
positive residuals of about 5°, 

(3) . The EURA-NOAM pole is shifted northward to 81°N, 11S°E, well 

outside the RM2 95% confidence ellipse. Consequently, the 

f 

variation in rates along this boundary does not match the 
• observations. 

(4) The AFRC-EURA pole is shifted westward to 12°S, 38°W. Such a 
pole implies compressive motion along the entire Azores- 

... Gibraltar Line. As noted above, this consequence is in direct 

" “ conflict with' the extension observed on the western portion 
-r ’ of this boundary. - ” ‘ * .“• 

We conclude that a non-zero NOAM-SOAM angular velocity is required 
by the revised data set. To derive RM2, we" adopted the convention of 
Paper 1 and partitioned the AFRC-NOAM and AFRC-SOAM data sets at-lS^N, 
where the distance” between the Mid-Atlantic Ridge and the West Indies 
Arc is least. ’ ’ ' ‘ ' ' ’ ■’ ■ 



Th ls-er-oupltl B ~af fords -afl'fexcsllyftt fit~€6"£hg~d/tfea'^Mong- ttie— " • 
-AFRe-rSOAlt'bocltidaty: ^Figure =16) i ^ne'daEutii^^on tMs^boutidaxy deserves 
particular mentlau; -’Elttrelm and‘'Eoing-(1975y~'Have:iinappe4“a"rccefttp" 
apparently continuous fault within the Veaia fault 'aohe, - Their data 
yield .a remarkably well determined azimuth of- relative motion; we 
assigned this datim an uncertainty of " ±2° j the loWeSf given to any -- --- 
direction datumr Its residual computed from RME-is-nrily 0.4°. In 
contrast, the residual -computed from the model x?ith h single Americas - 
plate is nearly -5°. 

Although some motion is required, the KOAM-SOAM angular velocity 
vector is not precisely constrained. This is indicated by the large 
confidence ellipse associated with , the pole (Figure 5) . It is also 
evidenced by the fact that the EM2 pole is nearly 30° north of ‘the RMl 
pole, completely reversing the sense of motion predicted along the 
boundary postulated to lie somewhere between 10°K and 20*^N. Discussion 
of the inferred relative motion may be found .in a later, .section. 

- • ■ - . ^ ^ . ’'z.. - . . " ' * ' ' 'I. r 

Caribbean plate motion. Although a Caribbean plate was not included 
in the RMl model derived in Paper 1, the topic of Car-ihbean plate 
motion was treated in detail by Jordan (1975) . He derived a NOAM-CARB 
angular velocity vector using -a spreading rate of 2.2 cra/yr' across the 
llid-Cayman Rise estimated from topographic decay (Holcombe et al. , 1973), 
For the present study, we were fortunate to have available a much more 
reliable rate (2.0 ± 0.4 cra/yr since 2.3 My B.F.) determined from a 
magnetic profile across the Mid-Cayman Rise by Macdonald and Holcombe (1978) 


_This rate is essentially identical to the previous estimate* Four slip 
vectors from the Molnar and Sykes (1969) set used by Jordan (1975) were 
deleted, one from the West Indies Arc, because it may lie south of the 
CAKB-NOAM-SOAH triple junction, and three from Hispaniola and the Puerto 
Hlco trench, where the data show internal scatter and the stress and 
strain fields are complex (Jordan, 1975) . A slip vector for, the 1976 
Guatemala earthquake (Kanamori and Stewart, 1977) was added. The changes 
to the direction data shifted the KOAM-CARB pole northwestward from the 
position computed by Jordan (1975). It can be seen from Figure 5, 
however., that this shift is in the direction least constrained by the 
data, as indicated by the orientation of the RM2 confidence ellipse. 
Jordan's pole lies within this confidence ellipse, and the difference 
between these poles is not resolvable by the present data set (Figure 17) 

The CA3XB-S0AM pole is also shifted with respect to Jordan's 
solution, but, again, the shift is along the major axis of the error 
ellipse* This pole is unconstrained by data along the CARB-SOAM 
boundary, so its 95% confidence ellipse is quite large. The change in 
its ~location-ref lent s. the shifts in both the NOAM-SOAM and NOAJI-CARB 
poles. Nevertheless, Jordan's conclusion that a component of north- 
south motion exists along this boundary is unaffected (Table 5), 

* • ^ ' 1 

" The Bouvet triple .junction, RMl did not predict correctly the 
relative motions of SOAM-ANTA and AFRC-ANTA (Forsyth, 1975; Sclater 
et al, . 1976). In Paper I, these boundaries were very poorly constrained 
by data, but this deficiency has been remedied by a number of recent 


special studies (Table 1). KM2 provides an excellent fit to the data 
around Bouvet triple junction (Table 1, Figures 16 aud 18) , whereas KMl 
performs miserably. Three explanations for this discrepancy were 
investigated: 

(1) RKL is located in a local minimuia of the fitting function 
manifold. This possibility can be dismissed; inverting the 
BKL data set with EM2 as a starting model yields the published 
BHl solution. 

(2) The SOAM-ANTA and AFRC-AHTA vectors are very sensitive to 
small errors in the BMl data set. This possibility can also 
be excluded; the epror ellipsoids for these vectors are 
actually quite small (Paper I, Table 2, Figures 5 and 7). The 
prediction error computed from the EMI variance matrix is 
much smaller than the RMl misfit to the new data. If the new 
data along the SOAM-ANTA and AFRC^’AHTA boundaries are exluded 
from the revised data set, a solution similar to RMl is obtained. 

(3) The global data set is inconsistent with the plate geometry 

assumed by KMl. 

•* 

Hypothesis (3) is our preferred explanation and was in fact advocated by 
Forsyth (1975) in his original study of this problem. For reasons detailed 
below, wc believe that the data sets for plate motions about the Indian 
triple junction are Inconsistent with our modpl, and we ascribe this 
Inconsistency to internal deformation within the Indian plate. 



25 


Plate motions in the Indian Ocean. This brings us to the major 
difficulty we encountered in constructing BU2: as pointed out by 

Jordan et al. (1976) and Minster and Jordan (1977) « each of the three 
legs of the Indian triple junction are populated by internally . 
consistent data, but the three best-fitting vectors sum to a vector 
(the closure vector) significantly different from zero (Table 3, Figure 6). 

The AFRC-AMTA boundary is densely populated by good observations. 

The 6 rates, 6 transform faults, and 11 slip vectors along this boundary 


constrain the angular velocity vector very well. The most important of 
these data is the well mapped Melville transform fault = 0.53) near the 
northeastern end of the boundary (Engel and Fisher, 1975), which controls 
the latitude of the pole. RM2 performs close to optimally along this 
boundary (Figure 18). 

As noted by McKenzie and Sclater (1971) , the transform faults along 
the Central Indian and Carlsherg ridges tightly constrain the INDI-AFRC 
pole, and these constraints have been strengthened by improved bathymetry 
(Engel and Fisher, 1975). As shown on Figure 19, there is a minor 
discrepancy between the rate data and the transform fault azimuths: the 

northernmost rates are too large by a few tenths of a cm/yr. In an 
effort to fit these rates, the best-fitting vector skews slightly with 
respect to the T.F. data, and EM2 is actually a better fit to the azimuths 
than the BFP. However, the Carlsherg Ridge is opening slowly and lies 
close to the magnetic equator; the magnetics along this boundary are 
not of exceptional quality (McKenzie and Sclater, 1971), and we are 
not disturbed by this slight misfit. 

+ A 

ft 
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The problem of data inconslsteni’y is evident along the Southeast 
Indian Ridge. The data are not quite as good along this boundary, but 
they determine a BFP and angular rate which constitute an acceptable fit 
(Figure 20) . KMl fits these data very well, but BM2 fits poorly; the 
RM2 pole is significantly different from the BFP (Figure 6) and does not 
match the gradient in the spreading rates. The situation is now clear: 

RMl satisfies the INDI-AFRC and IHDI-ANTA data, but misses badly along 
the AFRC-ANTA boundary; BM2 corrects the misfit, but then does not 
satisfy the INDI-ANTA data. The most comprehensive local study of this 
triple junction was published by McKenzie and Sclater (1971) . Their 
instantaneous motion model is also shown on Figures 18-20. It is 

different from either Kt-ll or RM2 but does not constitute a better solution. 

' ' * 

The motion of Arabia.. In the Gulf of Aden, the rates obtained by 
Laughton et nl. (1970, Table 1) are used directly. These data show very 
little scatter and are fitted by Bi'12 very well. The only other data used 
in the inversions are two rate estimates in the Red Sea (Allan and Morelli, 
1970), and these are also well fitted. Because* of the mediocre quality 
of the azimuth data, and the variety of the possible Interpretations 

of Red Sea tectonics (e.g, LePichon et al» , 1973), we did not attempt to 

^ * - 

model the northern Red Sea in this work. Since the Arabian plate is 
unconstrained along its other boundaries, the RM2 and best-fitting 
ARAB-AFRC vectors are identical. 
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The Indian Plate Problem 

Although KM2 is a very good fit to the data set as a whole, we have 
not been able to fit the Indian Ocean data satisfactorily by an 
RM2-type model. These discrepancies may simply result from bad data, 
contaminated by systematic observational errors we do not understand. 

Ve are aware that data bias is the probable explanation for the misfit 
to the Aleutian slip vectors; in Paper I, we attributed this misfit, 
evidently incorrectly, to .internal defo.rmation within the North American 
plate. The distance of systematic errors in the Indian Ocean data ob- 
viously cannot be ruled out at this time. However, because its impli- 
cations are important, an alternate hypothesis — internal deformation 

within the major plates — deserves investigation. 

» • 

In BM2, Indian Ocean tectonics are modelled by three plates, 

ANTA, AFRO and IHDI, There is no geological or seismic evidence for 
deformation within Antarctica; in fact, the intraplate seismicity of 
Antarctica appears to be the lowest of any major plate (e.g. Tarr, 1974). 
In contrast, both the African and Indian plates are characterized by high 
Intraplate seismicity, and observations of significant post-Miocene 
intraplate deformation have been' reported (e.g. McKenzie et al. , 1970; 
Sykes, 1970b; Eittreim and Ewing, 1972). 

To investigate hypothetical intraplate deformation, we have chopped 
these plates into two pieces and modelled each as a rigid entity, as we 
did for NOAM and SOAM. This procedure is obviously unsatisfactory for 
representing widely distributed strain, and we are impliclty asstjmlng 
that most of the deformation is localized within a relatively narrow zone. 


Deformation of the Afr!ican plate . Active extension across the 
African Rift valleys is well_ document ed _(a.g. McKenzie et al ., 1970; . 
Maasha.and Molnati 1972; LePichon et al. , 19731, ._Tq„ t§st the hypothesis 
that the_RM2 misfit along _the ,INI)I-AMXA_bpunda^ stems from ignoring 
this deformation, another global . inversion was .performed . The data 
along the African plate boundaries in . the Red Sea and west of 20°E.were 
assigned to a Nubian plate (NUBI), and the data east of 40°E were assigned 
to a Somalian plate (SOMA). We arbitrarily assumed that the position of 
the NUBI-SOMA-MTA triple junction is someswhere between 20®E and 40 “E, Since 
we did not feel justified in specifying its position more accurately, 
the 10 data along the Southwest Indian Ridge in this interval were deleted. 
As-. eJqjected, the resulting model is a better fit to the data set than RM2. 

In particular, the INDI-ANTA angtilar velot;ity vector is very Close to 
the best-fitting solution in Table 3, and the fit to data along this 
boundary is much improved. However, the I'esulting SOMA-NUBI pole is at 
43°S, 48°E and the angular rate is 0.17°/My, which implies east-west 
compressive motion across the African Rift valleys at a rate exceeding 
1 cm/yr! This prediction clearly contradicts the geophysical evidence. 

If a non-zero component of extension is imposed on this boundary, the fit 
to the INDI-ANTA data set is degraded with respect to RM2. 

Therefore, problems with RM2 in the Indian Ocean cannot be remedied 
by. simply postulating internal deformation in Africa, because the resulting 
model violates other constraints. Although the evid.ince for extension 
across the African Rift Zone is compelling, we have not been able to 
successfully resolve this motion in our global modelling studies, a 
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conclusion also stated in Paper I. 

In a recent parallel study. Chase (1978) has produced a global 
plate model which predicts opening of the Rift valleys. The differences 
between his model and the model described above are evidently due to 
differences in the Inverted data sets. We note that Chase's poles do 
not provide a satisfactory fit to our data set along the RM2 AFRC-ANTA- '' 
boundary. Also, the misfit to the INDI-ANTA data set described for 
RM2 is a feature of his solution as well. 

Deformation of the Indian plate. The hypothesis that the Indian 
plate is deforming is suggested by two aspects of the RH2 fit discussed 
in the previous pages: RM2's performance is unsatisfactory along both 

the IHDI-AMTA and INDI-PCFC boundaries. To test the hypothesis that 
INDI* deformation is responsible for these discrepancies, the western 
portion of the Indian plate (WIND) was separated from the eastern 
portion (AUST) . Six INDI-ANTA data within a transition zone between 
90°E and 130°E were deleted. Data on the Indian plate boundaries west 
of 90*^E were assigned to WIND and data east of 130° were assigned to 
ADST. With this configuration, the global data set was inverted. The 

resulting AUST-WIND angular velocity vector is labelled "A" in Table 4, 

* - • * 

Again, introduction of more model parameters permits'^a better fit to the 
observations: The remaining data along the Southeast Indian Ridge are 

satisfied, and the AUST-PCFC pole lies within 2® of the INDI-PCFC BFP of 


Table 3. 
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Prom Table 3 we can estimate the hypothetical AUST-inni) VECTOR inde- 
pendently of the data along the Southeast Indian Ridge. Deformation of 

the Indian plate can be approximately described by „the closure jyector 
of the circijtit TCElTO-APkc-AN’^-PCFC-AUST. 'This vector may -be calcinated 
using the best fitting angular velocity vector for each boundary traversed 
by the circuit. The result is net unique since the PCFC-AUST rate is not 
constrained, and a one parameter family of closure vectors is therefore 
generated. To specify a member of this family, we arbitrarily chose to 
minimize the relative velocity of AUST with respect to WIKD at a point 
along the Ninetyeast Ridge, Numerical experiments show that the result 
is quite insensitive to this point’s location. The derived angular 
velocity vector is labelled "B" in Table £■. 

In view of the uncertainties involved (and the ad-hoc criterion 

► * 

used to construct vector B), the two solutions in Table 4 are remarkably 
similar. Both imply slow compressive raot.ion between WIND and AUST in a 
NW-SE direction. 

Our modelling procedures do not require the existence of a specific 

boundary separating the Indian plate into two portions. However, we 

speculate that any deformation within the Indian plate may in fact be 

localized in the vicinity of the Ninetyeast Ridge. This linear feature behaved 
* * 

as an active transform fault in the Cretaceous (eg, McKenzie and Sclater, 

1971; Schlich, 1975; Sclater. et al, . 1976), and, although it has been 
commonly considered to be quiescent during recent times, Stein and Okal 
(1977) have suggested that it i!» now the site of significant seismic 
and tectonic activity. The nature of this tectonic activity is undoubtedly 
complex, but Stein and Okal argue that the bottom morphology and 
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eeiainiC;, source ffiechanisro? a re , consistent with IIW-SE t ompression .in the 

region, in agreement with the angular, velocity vectors in. Table 4. Vector 

A predicts a rate of deformation of about 1 cm/yr,- computed at IS^N, 

■"8 

90"E, This rate is equivalent to- a strain rate of 10 /yr, if. the. 
deformation were distributed over a sone 1000 km wide, and is grossly 
compatible with the level of regional seismicity (Stein- and Ohal, 1977) . 

It summary, the hypothesis that deformation is occurring within 
the Indian plate suffices to resolve, the. difficulties encountered in 
fitting the instantaneous relative motion data. Although the nature of 
this deformation remains speculative, at least a partial localization of 
the deformation in the vicinity of the Ninetyeast Ridge is suggested by 
other observations. We note that, if extension across the African Rift 

Zone is incorporated into the plate tectonic model, deformation w'ithin 

► • 

the Indian plate predicted by the model will be greater* 


: Along plate. boundaries-TAern data* aru' no travailablerorr where ‘inter- 
pretation is hindered by geological complications, EM2. provides a useful 
basis for predictions azuL comparisons of global motions with localr field 
evidence. Ue discuss here a few selected examples. In this* discussion, 
prediction errors were calculated using the bilinear form described- by. 

Jordan (1975). - . * - 

. Central California. Because of possible- bias associated with ex- 
tension in the Basin and-Range Province, data along the San Andreas - fault 
system were not used in the inversion (Figure l).*In central California 
EM2 predicts a rate of relative motion between the Pacific and North Ameri- 
can plates of 5.6 ±0.3 cm/yr (Table 5). Based on geological evidence. Hall 

# * 

and Sieh (1977) estimate a slip rate of 3.7 ± 0.3 cm/yr along the San 
Andreas in central California, averaged over three tnillenia, which is 
identical to Thatcher’s (1977) geodetical estimate of 3.7 ± 0,2 cm/yr. 

Ci i'li-';’ evidence irelios a sir.rllai* rntc over the past 10 My (eg, 

Huffman, 1972). This comparison suggests that a significant fraction of 
the PCFC-NOAM motion is taken up elsex^here. Some of it may possibly be 
accommodated on fault systems west of the San Andreas. For example, 

Weber and Lajoie (1977) conclude that right- lateral slip has occurred 
along the San Gregorio fault none during tlie last 200,000 years, with 
a rate ranging from 0.6 to 1.3 cm/yr. The discrepancy between the 
observed and predicted rates may also be attributed to deformation 
distributed \d.thln the Basin and Range Province. Thompson and Burke (1973) 




estimate that the Basin and Range tmderwent 100 km of extension In H55W 
direction during the last 15 My, equivalent to an average of 0.7 cm/yr 
of right-lateral motion in a direction parallel to the San Andreas. 

The comparison of observed and predicted azimuths also suggests' 
active deformation within the western U.S,: Between the Carrlzo Plain 

and Hollister, the San Andreas fault exhibits a well defined azimuth of 
N41"W ±2°, whereas the direction of relative motion calculated from 
BM2 is U35“U ±2“ (Table. 5). These two values can be reconciled by 
postulating about 0.8 cm/yr of EtJ extension between central California 
and the stable North -American platform to the east. Thompson and 
Burke’s (1973) model implies an average rate for EW Basin and Range 
extension of 0.5 cm/yr. Furthermore, Clark and Lajoie (1974) estimate 
a horizontal displacement rate' of 0.7 cm/yr along the Garlo’cfc Fai:lt 
during Holocene time. Such agreement may be fortuitous, but we consider 
it to be support for Davis and Burchfiel’s (1973) suggestion that the 
Gnrlock Fault is a major Intracontinontal transform structure. 

Relative motion of North and South America . As argued above, 
relative motion beti^een North and South America is required by our data 
set. Figure 5 and Table 2 indicate that the NOAM-SOAM vector is poorly 
constrained and a wide range of possible relative velocities are allowed 
hy the data. Very little direct evidence for this relative motion exists 
and the movement could be distributed across a broad zone between, say, 
10*N and 20**N. Since the relative velocities are predicted to be small, 
the deformation may be largely aseismlc* However,' some seismicity does 



exist. For example, a magaitude 6*2 earthquake occurred October 23, 1964, . 

at W.S^N, 56*1®W, The mechanism for this event is consistent with right- 

lateral strike-slip motion in a direction N5S“W (Molnar and Sykes, 1969; 

J* Sorel, 1975, personal communication), whith does not disagree with 

the KM2 prediction of Nyi^W ± 58“ (Table 5). It is, however, inconsistent 

with the Rttl model, which predicts left-lateral motion. 

Southern Boundary of the Caribbean plate . EH2 predicts a component 

of NS convergence across the CARB-SOAM boundary. Although the rates are 

« 

somewhat higher, the azimuths for CAEB-SOAM motion are almost identical 
to those deduced by Jordan (1975) using the EMI model. Consequently, 

Jordan's conclusions concerning motions along this boundary are substantiated 
by this study. They are also supported by Ladd’s (1976) model of tertiary 
plate motions. Direct evidence for NS compressive motion has been obtained 
by Taiwan! al . (1976) from an analysis of multichannel seismic reflection 
records from the south margin of the Venezuelan Basin and by Rial (1978) 
from a study of local mechanisms in Columbia and Venezuela. No such com- 
pression is predicted by a model which assumes a single American plate. 

We take this to be an additional argument in favor of modelling NOAM and 
SOAH as two separate plates with a zone of decoupling between 10“N and 20“N. 

... Jordan's (1975) portrayal of the tectonic relationships in the 
Panama Basin is also compatible with EM2. The EM2 ' COCO-NAZC pole lies 
north of the equator, and the Panama T.F. , as mapped by Lonsdale and 
Klitgord (1978) , closely approximates a small circle about this pole', 
even though it was not used in the inversion. Thus, EM2 is consistent 
with the hypothesis that the Panama Basin east of this transform is not 
acting as a separate plate, as suggested by Jfolnar and Sykes (1969) and 







Idnsdale kiitgOfTicl9^^ pa^'Vf' the liizca'piaia^ 
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asdinuths of relative motion, are nearly identic^ (Table 5) and are con- 
sistent with the hypothesis that the motion is accoimnodated by a left- 
lateral transfoinn fault along the southern continentiL margin of Panama 
(Jordan, 1975), • 

Subductxon ~ia~gb'ntHern “Cfille . Seismic activity along the Chile 
trench decreases sharply south of the NAZC-A3ST&-S0AM triple junction 
'(Tarr, 197*4). Pew earthquakes (only one with nij^ > 6) have been 
reported in this region between 1963 and 1975. The predicted convergence 
rate between ANTA and SOAM is only 2.1 ± 0.2 cm/yr. (Table 5), 6.7 cm/yr 
less than the subduction velocity north of the triple junction and 30% 
lower than the EMI prediction. Yet other convergence zones with com- 
parable rates such as the West Indies Arc or the South Sandwich Trench 
are significantly more seismically active. If our model is 
correct, then subduction in Southern Chile takes place largely aseismically , 
or this boundary constitutes an extensive seisHiic gap. 


The Owen Fracture Zone . The Owen Fracture Zone represents the-INDl- 
AEAB boundary (e.g., McKenzie and Sclater, 1971) and exhibits only weak 
seismicity. As shown in Table 5, RM2 does predict a low rate of relative 
motion between these two plates, but the predicted azimuths do not agree 
'Well with the observations. At 14 ®N, Laughton’s (1970) bathymetric map 
indicates an azimuth of N30®E for the Owen fracture zone, compared with 
the model value of N55®E ± 14“ , and at 22“K, a fault plane solution by 
Sykes' (1967) has a slip vector orientation of N50“E, versus a model value 



^eaqjaata. sHgge st^th at the 
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tli e^JLiiye rsion with .IJSPI ^ ej)arat^ int.Q_WIl^^^^^ST^^eacrib€t4 ahoy e ,_^. , 
yields an WIH^AIl&B pole positioned 3°^ north o£ the: JtM2_ poJLe. Tha_ - 
azimutV calctd^^ M°N». 59°E JLs ^^t.ter.jasre^ent .with, the^ 

observations j e.1 though the azimuth calculated at 22 *N, 62*E is nearly 
identical to tliat for EH2. _ . 
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.*r£ rf '^r. £ * - V.=ri-=- Absolute Motions . - ■ , ■.• - . 

The BM2 geohedron (Table 2, Figure 2) compietely describes the 

relative motion model. To specify an 'absolute* reference frame, we 

need only to choose an origin in angular velocity space. A particular 

frame of interest in discussions of plate dynamics is one fixed with 

respect to the average position of the deep mantle, assumed to be rigid . 

« 

or at least to have typical internal motions much slower than the motions 
of the plates; we refer to this frame as the mean mesospheric frame. 

In Paper I we constructed an absolute motion model (AMI) based on 
the Wilson-Morgan fixed hotspot hypothesis and concluded that this hypo- 
thesis was consistent with the available Instantaneous motion data. 

However, we noted the difficulties in estimating rates and directions of 
hotspot migration that are compatible with the short time intervals 
appropriate to the relative motion model, especially for hotspot traces 
on the slower plates; Because of these difficulties, we are intrinsically 
limited in our ability to construct more refined tests of the Wilson- 
Horgan hypothesis and to discriminate among various instantaneous absolute 
motion models using hotspot data. 

To investigate this limitation, we have derived an absolute motion model 
by again inverting hotspot data, but restricting the data set to include 
only those constraints on hotspot migration pertinent to the last 10 Hy. 

This time span is really the minimum interval for which good hotspot data can 

be obtained, although it exceeds by over a factor of three the mean 
averaging interval for the relative motion data. The azimuths of nine 
hotspot traces and the rates for five were chosen on the basis of this 
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criterion (Table 6) . The data set is dominated by the information from ■ 
Pacific island chains; no Atlantic or Indian Ocean hotspots were employed. 
The rate at Hawaii represents our interpretation of the K-Ar ages between 
Hawaii and French Frigate Shoals summarized by Dalrymple ^ (1974) . 

For four other Pacific archipelagos the K-Ar ages of Duncan ^ a^. (1974) 
and Duncan and McDougall (1974> 3976) have been used. Azimuth estimates 
for the traces were obtained from bathymetric charts , and -the rate esti- 
mates were projected along these directions. The mean rate estimates for 
individual island chains have formal standard errors of about + 1 •cm/yr , 
(Duncan and McDougall, 1976) , but these have been increased to allow for 
possible errors due to biased sampling. (We note that, since vulcanism 
may persist at a given site for millions of years, a systematic failure 

* “u 

to sample the oldest rocks generally results in rates biased to high 

values.) The other data in Table 6, hotspot azimuths from the COCO, NA2C 

and NOAM plates, have been taken from Paper I. 

The data- set in Table 6 was inverted to obtain an absolute motion 
model designated AMl-2 (Table 7, Figure 2). In the inversion the relative 
plate velocity vectors were fixed at their EM2 values, but the uncertain- 
ties in the RM2 model, represented by its variance matrix, were incorporated 
into the calculation of the variance matrix for AMl-2. The model Is a 
very good fit to the selected data set: only one datum has a residual 

exceeding its assigned error (the azimuth of the Marquesas) , and the rate 

data are all fitted to within 1 cm/yr. Thus, the results of this experi- 
ment give us no cause to challenge the Wilson-Morgan hypothesis. 

But, even supposing the Wilson-Morgan hypothesis is valid, which we 
have not proved, with what precision can the motions- of the plates in the 
mean mesospheric frame be predicted by the hotspot data? The answer to 
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this, question is indicated by the standard errors of estimation listed in 
Table 7, Although’ the absolute velocities of the fast--nioving oceanic 
plates {e.g. PGFC) have relative errors which are small, the relative errors 
for the slowly-moving continental plates (e.g. EURA) are quite large and in 
some cases exceed 100%. Hence, the absolute motion directions of several 
plates, particularly AHTA and EURA, are not usefully constrained by the 
hotspot data used in this experiment. For example, at the position of . 
Iceland the motion of EURA with respect to the mean mesospheric frame is 
predicted by AMl-2 to be N83°W at 0.4 cm/yr, nearly diametrically opposed 
to the direction of the Wyville-Thompson Ridge, the presumed hotspot trace. 
Sut no significance should be assigned to this discrepancy, since the formal 
ipxediction errors (la) are + 1'62° and + 0,8 cm/yr, respectively, and since 
. . the actual azimuth of the Iceland hotspot trace over the last 10 My is not 

really l^nown (Paper I, p.566). 

With these large uncertainties in mind, it is interesting to compare 
the hotspot model with absolute motion models based on other criteria. 

Three such alternate models are listed in Table 8 (see also Figure 2) . 

AMO-2 is the unique absolute motion model constructed by requiring that 
the lithosphere as a whole possess no net rotation, a criterion discussed 
and applied in Paper I and by Lliboutry (1974) and Solomon and Sleep (1974). 
AM2-2 corresponds to Burke and Wilson’s (1972) hypothesis that the African 
plate is stationary with respect to the mantle, a criterion endorsed by 
Duncan and McDougall (1976) on the basis of Pacific hotspot data, AM3-2 
.conforms to Jordan’s (1975) suggestion that the Caribbean plate Is fixed 
in the mean mesospheric frame, pinned in position by its two bounding 
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subdue tlon zones. 

The predictions of these absolute motion models are compared with 
the selected hotspot data in Table 7. The Pacific poles for all of the 
absolute motion models are similar (Table 8) , and the azimuths of the 
Pacific island chains are essentially equ ally well fitted by each. However, 
the Pacific rate data and the azimuth data from the other plates do 
provide some discriminants. AMO-2 appears to be inconsistent with the 
rate data; its values are significantly less than those observed. AM2-2 
is a good fit to the Pacific data, but it is a poor fit to the azimuth 
data for the other three plates. AM2-3 provides a good fit to the azimuth 
data, but its Pacific rate is slightly low. 

The alternative absolute motion models can be compared directly with 
AMl-2 in model space using the computed estimation errors. , Let m^ be the 
model vector representing AMl-2 and let in‘ be any altenative absolute motion 
model. .Define the quadratic form 

F » (m - m')^ ♦ V 

I 

2 

where V is the complete variance matrix for m* Then, if F > (1.96) , 

^ lies outside the AMl-2 95% confidence hyperellipsoid, and one can accept 
the conclusion that the expected value of m (of which m is only an estimate) 
is different from m’ at the 5% risk level, (Of course, this statement assumes 
that normal statistics and our linear approximations are applicable and that 
V is known exactly, which is not strictly true; it nevertheless provides a 
workable basis for making statistical decisions.) For models AMO-2 and AM2--2, 

F equals 12.4 and 10.9, respectively; we conclude that these frames are 
significantly different from the hotspot frame. For AM3-2, F equals only 3.1, so 


the hypothesis . that the Caribbean plate is fixed in the hotspot frame 
zzr...- ■ 

cannot^be rejected. We. note that the frames corresponding to ANTA fixed 
(F = 0.5) and EURA fixed (F = 0.8) are Indistinguishable from the hotspot 


frame as well. 


It is also interesting to compare AM, the absolute motion model 
derived in Paper I, with AM-2. Both models were obtained by the inversion 
of hotspot data, but, in the case of AMI, no rate data were used and a 
much larger, more globally distributed set of hotspot azimuths were fitted. 
As a consequence, the averaging intervals for the AMI data are generally 
greater than 10 My and more variable. Although the AM and AM-2 Pacific 
poles are similar, the AMI rotation rate (0.83 °/My) is less than that of 
.^iM-2 (0.97 °/My). For AM, P = 339. This very large value is indicative 
©f the fact that RMl and RM2 are significantly different relative motion 

, t 

models, in that RMl lies well outside RM2’s 95% confidence hyperellipsoid. 

A model derived by adding to SM2 the AM PCFC absolute rotation vector 
yields F‘ = 10.0 and is inconsistent with the data set in Table 6. 

The resolution of absolute motions by the hotspot data is obviously 
degraded if the possibility of a non-rigid hotspot geometry is allowed. 
Several authors have concluded that, averaged over geologically long 
periods of time (> AO My) , hotspots have relative velocities with mag- 
nitudes on the order of 1 cm/yr (Morgan, 1972; Burke et al. , 1973; 

Molnar and Atwater, 1973; Molnar and Francheteau, 1975). In some sense, 
oiir 'conservative assigr-.ent of large errors to the hotspot data in Table 6 
may account for the uncertainties generated by small random motions among 
the. hotspots, but appropriate caution in Interpreting any hotspot model 


must be exercised until better data and more rigorous tests are available. 



Nevertheless, several previously published conclusions regarding 
present-day absolute motions appear to be warranted; ■ the'se are common 
to all of the models in Table 8; 

(1) Plate speeds correlate negatively with total continental 
area (Paper I) . 


(2) Plate speeds correlate positively with the fraction of plate 
boundary being’ subducted (Jordan and Minster, 1974; Forsyth 
andTJyeda, 1975). 

(3) Plate speeds correlate positively with geographic co-latitude 
(Solomon et al. , 1975) , 

Simple mechanical models have been formulated to exp3ain the first two of 
these correlations (Forsyth and Uyeda, 1975; Solomon et al. , 1975; 

Kaula, 1975), but their true dynamical significance is still quite 
speculative. For example, SolOmon et al. (1977) have sugges*ted that 
these aspects may have very little to do with dynamics; they argue that 
the absolute plate motions characteristic of Tertiary time exhibit 
none of the correlations stated above. Although we eye their reconstruc- 
tions and modelling assumptions with some skepticism (cf . Jurdy, 1977) , 
we agree that more refined tests of the mechanical models must be 


formulated. 






ORIGINAL PAGE’ IS 

■ , OF POOR QUALITY 

Perspective 

BM2 is a significantly better representation of present-day plate 
motions. thau. EMI. In a recent parallel study. Chase (1978) has presented 
a global plate motion model generally quite similar to RM2. Some significant 
differences between these two models do exist — most ascribable to differences 
in data selection and interpretation— but the overall agreement is encouraging. 
These studies should be viewed as ever more rigorous tests of the plate tec- 
tonic hypothesis. We continue to be impressed by how well the large data 
sets (330 members in Table 1) are described by simp3c models with very 
few parameters (30 for EM2). 

We have noted, however, several problem areas where the plate model 
does not adequately fit the observations. These discrepancies deserve 
special scrutiny; they may be the manifestations of tectonic processes 

* f 

or other physical phenomena not now understood. For example, if our 
hypothesis that the Indian plate is not behaving rigidly is confirmed by 
better data in the Indian Ocean, then several questions must be addressed. 

How is the deformation distributed within the plate? What is the nature 
of the forces driving the deformation? Consider the hypothesis that the 
deformation is localized in the vicinity of the Kinetyeast Eidge: then 

a situation exists where, on two opposing plates at approximately equal 
distances from their common boundary (a spreading center) , there are two 
HS-trending zones of deformation, one extensional (the African Rift) and 
one compressional (the Ninetyeast Ridge) . This unusual configuration 
should provide a strong discriminant for force-balance models of the 
sort proposed by Forsyth and Uyeda (1975) , Solomon et al. (1975) and 
Richardson et al. (1976). Of course, more data are required before this 



hypothetical situation can be accepted as .reality. 

Throughout .the bulk -of this- 

tectonics haye been carefuyy. avoided. It- is clear_ that, in most jr.egipns 
of intracontinental .defonn.atipn, the pla-te ®odel^_has only limited utility, . 

However, global plate motions do provide the displacement boundary 

conditions -required to understand the .kinematics and dynamics of tectonics 
In complex reigions (e.g. Molnar and Taponnier,. 1975). These complex 
regions include not only the continental Interiors, but also .zones of 
deformation along the continental margins (e.g. Jordan, 1975) and even 
boundaries between the oceanic plates themselves. It is possibly 
complexities of this latter type which- are responsible for the difficulties 
we experienced in obtaining closure about the Azores triple junction. 

Unlike the relative motions, the absolute motions of plates in the 

* 

9 

mean mesospheric frame cannot be precisely constrained. Absolute motion 
models have been derived from a number of kineinatical hypotheses, and, 
although' they are grossly similar, significant differences among them do 
exist'. In our opinion, model AMl-2, with its attendant uncertainties 
(Table 7) , represents the most satisfactory description available from 
the present obseirvatlons. Based on these absolute motions, a number of 
empirical correlations appear to be warranted, but how these correlations 
relate to the fundamental forces driving the'plates is only speculative. 
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Appends x . 


-In—tha-intetpretation of- earthquake mechanisms along stibductlon 

boundaries, most authors assume that the direction of relative plata 

motion is given'by the horizontal projection of the slip' vector—(«-.g. 

Paper I). If the convergence is oblique to the trench axis, this 

procedure yields a biased estimate of the direction of relative motion. 

* 

Instead, the slip vector should be rotated into the horizontal plane, 
vhich requires correcting the slip vector azimuth by an amount- a- given by 


'cot(T - T ) 

o - orccot I p -- 1 + Tp - 


where T^., Pj,,and T^, are the azimuth and plunge of the poles of the 
fault plane and auxiliary plane, respectively. 

9 « 

Ihis correction was applied to the data from the Aleutian-Kuril, 
South American and Tonga-Kermadec Trenches, The statistical information 
is summarized below: 



1 “! “ 

NOMl-PCFC 

0.6" 0.3 

NAEC-SOAM 

0.9" 0" 

PCFC-INDI 

1.1" -0.9 

correction 

is clearly minor. 


I a max I 
2 “ 

2 “ 

4 “ 


omitting this correction does not give rise to a significant systematic 


bias in the data. 
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*For rates (MA) , units are cni/yr; for transform fault (TF) 
and .slip vector (SV) azimuths, units are decrees. 
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2. RM2 

Caohedron 









BeleCtve Rotetlon Veccor* 


error Rlllpaa** 

iBportonca DiaCrlbutloQ 

ri«M Pair 

Q 

®6 

♦ 

% 

u 



a 

ux 

®»la 






(*H) 


(°8> 

(dea) 

(®/My) 

< /My) 

(dcg> 

<deg) 

(deg) 

KA 

TT 

8V 

Total 

NOAM-PCFC 

48,77 

1.10 

-73.91 

1-94 

0.852 

0.025 

S71S 

1,30 

1.08 

0,405 

0.398 

0.694 

1.497 

CaCO-PCFC 

16,72 

0.89 

-107.39 

l.Ol 

2.208 

0.070 

S37E 

1.00 

0.63 

0,977 

0.272 

0.009 

1.258 

RAZC-PCFC 

56.64 

1.S9 

-87.88 

1.81 

1.539 

0.029 

8098 

1.91 

0.96 

0.669 

0.341 

0.038 

1.228 

£UtJl*PCFC 

60.64 

1.04 

-78.92 

'3.04 

0.977 

0.027 

5782 

1.31 

1.02 

0 

0 

0 

0 

nroi-PCFC 

60. 71 

0.77 

-5.79 

“ _1.83 

1.246 

- 0,023 

• S82S 

0.90 

0.76 

0 

0 

0.246 

0-246 

AHTA-FCFC 

64,67 

0.90 

-80,23 

2,32 

0,964 

0,014 

NS2E 

1,11 

0.75 

1,200 

0.811 

0,039 

2-050 

COCO-NOAM 

29. 8Q 

1.06 

-121.28 

■ 2,07 

1.489 

0.070 

S75S 
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0.99 

0 
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0.165 

AF&C-KQA.'t 

80.43 

1*57 

56.36 

35.29 

-D.2S8 

0.019 

KB6E 

5.88 

1.51 

0,851 

0.246 

0.091 

1,180 

EUEA-KDAH 

65,fl5 

6.17 

132.44 

5.06 

0.231 

0.015 

SliE 

6,36 

1.39 

1,055 

0.626 

0.366 

2,047 

HOAH'CARa 

-33.83 

9.19 

-70,48 

:2.76 

0,219 

0.052 

SUE 

9.42 

0.97 

0.952 

1-741 

0.253 

2,946 

COCO- CARS 

23.60 

1.48 

-115.55 

2.26 

1,543 

0.QS4 

5632 

2.24 

1.21 

0 

0 

0,1U 

0.111 
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47.30 

S-37 

-97.57 

4.57 
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0,056 

S19E 
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2.67 

0 

0 

0 

0 

COCO-KAZC 

S.63 

1.40 
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2.61 

0.972 
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1.40 

1.829 

0.732 
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25.57 
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-53,82 

6.22 
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7.22 
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0 

0 

0 

0 
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73,51 

11.75 

60.84 

48.86 
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0.038 
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0 

0 

0 

0 
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59.08 
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0.835 
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0,464 

0,464 
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0.608 

0.283 

1.058 

IHDI-AFRC 

17.27 
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3.44 
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0-024 

S65E 

2.72 
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0 
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0.222 
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0.009 

5422 

4.93 
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0.697 
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2.135 
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18.67 

1.16 

32.74 

1.41 
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0.011 

5622 

1.39 

1.10 

1.012 

0,135 

0.025 

1-172 










Totals 

14,273 

11,134 

4.593 

30,000 

plate 

n^sed Doves counter clockwise with 

respect 

to the second. Oncercalntlca are the 

standard 

devlatiens of 


Docginal dlstrlbutletia* 












**One efsma error ellipses are 

specified by the 

aslauth 

of the major axis 

lengths of the 

axes are geocentric angles. 
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Table 3, Beat-fitting 

angular 

velocity 

vectors for Individual plate boundaries ,* 





Relative Rotation Vector 


Error Ellipse 
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( /Hy) 

(deg) 

(deg) 

(deg) 

MA 

TF 
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Total 

KOAH-PCFC 

49.02 

1.34 

-76.05 

2.62 

0.865 

0.071 

K66E 

1,85 

1.34 

1.000 

0.615 

1,385 

3.000 

COCO“PCFC 

37.68 

3.39 -107.74 

1.84 

2,298 

0.317 

SlOE 

3.44 

1.35 

2.002 

0,965 

0.033 

3-000 

KAEC-PCFC 

55.64 

6.72 

-85.76 

5.40 

I.S27 

0.062 

N22E 

7.20 

1.62 

1.978 

0,928 

0.094 

3.0C0 

IKDI-PCFC 

55.71 

2.75 

-5,00 

3.42 

_ 

— 

519£ 

2.64 

1.79 

0 

0 

2.000 

2,000 

-T’ri r 


1.55 

-79.26 

3.40 

0.976 

0,022 

H41E 

1,92 

0.85 

1.600 

1,150 

0.050 

3. COO 


• \ * ■ 

U..51 

-11. 

60.61 

0.275 

r.n f-. 

N7 3E 

17.91 

3.38 

1,706 

0,977 

0,317 

3.0C0 


t , ■ * 

9,72 

i:7.37 

9.83 

0.252 


r : 

! ‘ . r L 

1.1) 

l.T 

r.772 

0.4-7 


Ka\H-CAK5 -34. lb 

9.28 

-70.40 

2,80 

0,225 


bl >- 


Lr,'# 7 

l.t 

1 Y s 

0 . . ' 5 

- ' 

COCO-HAZe 

4.90 

1.54 

-123.65 

2,71 

0.993 

0.071 

S69E 

2.70 

1,54 

2.003 

Cr.903 

O.U.^4 

3 , ^ 

APRC-SOAH 

62.98 

3.46 

-39.14 

2.46 

0,357 

0,010 

S06E 

3,48 

1.07 

1,407 

1.491 

0.102 

3.000 

IKDI-AFRC 

19.63 

1.7B 

44.12 

1.81 

0.612 

0,035 

S44E 

2,34 . 

0.76 

1*662 

1,318 

0 

3,000 

ARAB-AFRC 

30.82 

3.44 

6.43 

11.48 

0.260 

0,047 

S79E 

10,02 

2.93 

1. 989 

0.934 

0.077 

3,000 

AFRC-EURA 

35.71 

4.07 

-21.04 

0.95 

— 


SOIE 

4,07 

O.BS 

0 

0.793 

1.207 

2,OCO 

AFKC-AVTA 

7.93 

4,98 

-38.72 

4.77 

0,146 

0.011 

S43E 

6,70 

l.SI 

1,038 

1.666 

0.296 

3.000 

IKDZ-AirrA 

11.85 

2.76 

34.74 

2.99 

0,672 

0.013 

H6SE 

2,98 

2.71 

2, COB 

0.8)9 

0.153 

3,000 


*Sv«bQli> And conveocloiiA tho iAt&a *b Table 2. 
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Tablf> 6. A ccTDpArlaon becveen hotapoc dau <nd abaolutu motion modalf . 


Hot*poC 

Ut. 

Lgnfi* 

rt) 

Tract 

ri/tte 

Observed 

Azlnoth Race 

(deg) (co/yr) 

AKl-2 

ArlButh Rate 

<dcg) (co/yr) 

Importancfi Distribution 
Asinuths Racea 

AMO-1 

Azimuth 

(deg) 

Rate 

(cm/rr> 

AM2-2, 

Azliuth Rata 

(deg> 'Xeti/yr) 

AH3-2 

AaljKich 

(deg) 

Rate 

Cca/jn 

BawAtl 

20 

-X55 

HavalXao Zalanda 

pcrc 

H64W + 10 

10.0 ± 2,0 

H60H 

9.7 

0.292 

0,141 

K61V 


7.7 

K57U 

10.9 

K62U 

1.0 

Karqu^gaa 

*U 

-X38 

Karqueaao lalanda 

PCFC 

W43W + 15 

9.8 + 2.0 

W67H 

10.7 

0.078 

0.200 

H65U 

I 

B.l 

HfiZU 

11.6 

,I 

H67V 

9,5 

lAhltl-Hab£tU 

-la 

-148 

Society lalaadff 

vcyc 

W65W + X5 

11.0 + 2,0 

K64W 

X0.7 

0.084 

0.225 

H63W 


8,0 

K60W 

1: 11. S 

tf65U 

9.4 

Hacdonald 

-29 

-lAO 

AuatraX talanda 

rcFC 

H55W ± IS 

10.5 + 2*0 

H66H 

^ 10.5 

0.074 

0.239 

B63W 


7.6 

H60H f: 

|v u.a 
1 

H6SW 

9,1 

PiccAlrn 

~2i 

-130 

Plccalm-Cuablcr 

PL^C 

H65U + 15 

11.0 + 3.0 . 

K70H 

10.7 

. 0.065 

0.099 

H66V 


6.1 

H64W : 

11.3 

H66U 

9,4 

JuAQ de fucm 

. 46 

-X3G 

Cobb SeaoouoCa 


US4W + 15 

“ 

N47W 

5.5 

0.389 

- 

H49W 


5.0 

H43U 

: 7.0 

j 

H50U 

5,4 

CalApa^oft 

-1 

-92 

Cocoa Xldge 

C ; 

H45E + 1C 

- 

H46E 

8.5 

0.174 

- 

H52E 

1 ■ 

10.6 

K3Se 

r' 9,2 

1 ; 

H49£ 

9.5 

C«tap£goa 

-1 

-92 

Carnfegla Hldge 

yrc 

sasE + 10 

- 

H88E 

4.9 

0.480 

- 

HS3B 

» * 

7.2 

M71E 

'; 4.T 

K652 

S.f 

Tetlovatona 

45 

-110 

SnaVe River Plain 

r^A'i 

S60W ± 20 

- 

556W 

2,4 

0,460 

«a 

SSStf 


2.i 

H87W 

1 2.3 

'1 

S58tl 

2.4 
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T«fale 7. Modal AKl-2 


AbsoluCe RoCAclon Vectoi: Error Elllpvo 


riBCe 

0 

^0 
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ta 


^BS3t 




en) 

(deg) 

C^E) 

(deg) 

(®/Hy) 

( /My) 

(deg) 

(deg) 

(deg) 

AFKC 

IB, 76 

33.95 

-21.76 

42.20 

0.139 

0.055 

S73E 

40.43 

33.24 

ANtA 

21,85 

91-31 

75,55 

63,20 

0.054 

0,091 

H12E 

93.01 

56,12 

ARAB 

27.25 

12,40 

-3.94 

18,22 

0,386 

0.067 

S7tE 

16.38 

12,11 

CAM 

-42o80 

39,20 

66.73 

40.98 

0.129 

0.104 

N30E 

43.21 

23,90 

COCO 

21.89 

3.08 

-115.71 

2. SI 

1.422 

0.119 

532E 

3.35 

2.25 

EDRA 

0,70 

124.35 

-23.19 

146.67 

0.03B 

0.057 

S67E 

151-10 

116,90 

IKDl 

19,25 

6.96 

35.64 

6.57 

0.716 

0,076 

S25E 

7.16 

5.97 

KA2C 

67.99 

9.36 

-93.81 

9.14 

0.585 

0.097 

so:e 

9.37 

S.43 

KOAK 

-5B.31 

16.21 

-40.67 

39,62 

0,247 

0,060 

S57E 

23,12 

12.14 

PC3tC 

-61.66 

5.11 

97. 19 

7,71 

0,967 

0.08S 

SI6E 

5,23 

3.50 

50AH 

, -82.28 

19.27 

75.67 

85. 68 

0.285 

0.084 

H03E 

19.28 

11,38 


Table 8* Absolute isotlon ciodcis* 

?«cl£lc Rocotloo Vector 
Ifodel Klneaaclcsl Condition Long. Rate 

<“h) <“e) (“;Ky) 


AHO-2 

Ho net Totaclen 

-62.93 

111.50 

0-736 

AHl-2 

Best fit tc batapot dsta 

— 61. 6f 

97,19 

0,967 

AH2-2 

African plate fixed 

-59,15 

109.60 

1-043 

AM3-2 

Caribbean plate fixed 

-63,52 

104,4$ 

0.853 
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Figure Captions 


Figure 1. 


Figure 2. 


Figure 3. 


l-'iguro 4. 


Figure 5. 


Plate geometry and geographical distribution of the data 
used In producing model PM2. Circles Indicate seafloor- 
spreading rates, squares represent transform faults, and 
triangles slip vectors. Seven EURA-NOAJI data at high 
latitudes are not shovm on the figure. 

RM2 geohedron (stereo pair) . The geohedron depicts relative 
motions in angular velocity space (McKenzie and Parker, 197A) . 
Individual plates correspond to vertices. The z axis coincides 
with the rotation axis of the earth, the x axis is along the 
Greenwich meridian. Vectors representing the three reference 
axes have a magnitude of 0.3"/My. Open circle is coordinate 

« 4 

origin for AMO-2* Closed circle is coordinate origin for AMl-2* 
Histograms of normalized residuals for each data type, with 
sample size, sample mean and sample variance. The theoretically 
ideal Gaussian distributions with zero mean and unit variance are 
shown for comparison. Shaded area in lower histograms represents 
residuals for Aleutian and Kuril slip vectors, which show negative 
hiciy^, 

Poles for model RM2, with choir 95 percent (2a) coniidcnce ollipt^cs 
RMl poles and best fitting polos whore available (BIT, Table 3) 
are also shown* 

See Figure 4* 


Figure 6* See Figure 5 


Figure 7. In Figures 7 ~ 20, data and models are shown as residuals 
with respect to the predicted values calculated from the 
best fitting angular velocity vectors. Azimuths are measured 
In degrees in a counterclockwise direction, rates are in 
centimeters per year. Data symbols are the same as in 
Figure 1. Error bars are the subjective error bars listed in 
Table 1. The solid lines represent model RM2 (this study) 
and the dashed lines represent model EMI (Paper I) . Here 
the dashed-dotted line corresponds to the Bering-Pacific 
pole determined in Paper I. 

Figure 8. See Figure 7. 

Figure 9. See Figure 7. 

Figure 10. See Figure 7, 

Figure 11. See Figure 7. 
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Figure 12. See Figure 7. UK POUR QUALITY 

Figure 13. See Figure 7. 

Figure 14, See Figure 7. 

Figure 15. See Figure 7. 

Fifi;r<? 16. See Figure 7. 
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Figure 19. See Figure 18. 

Figure 20. See Figure 18. 
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